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1.  ABSTRACT 


The  heat  transfer  to  a  hemispherical  nose  cone  sutjec- 
ted  to  hypersonic  flow  conditions  is  calculated  for  three  Revrcii 
Numbers.  The  theory  is  compared  with  the  results  of  test  runs 
at  these  Reynolds  Numbers  on  a  stainless  steel  hemisphere  in  the 
Polytechnic  hypersonic  tunnel  facility.  Except  for  the  region  of 
highest  (theoretical)  heat  input,  the  comparison  is  good. 

Using  the  theoretical  heat  transfer  results  in  conjunc¬ 
tion  with  the  one-dimensional  transient  heat  conduction  solution 
through  a  spherical  shell  element,  temperature  histories  are  ob¬ 
tained  at  selected  points  in  the  interior  of  the  hemispherical 
nose  cone.  A  comparison  with  experimental  data  shows  a  satisfac¬ 
tory  agreement. 

It  is  concluded  that  the  theoretical  heat  transfer  va¬ 
lues  used  in  conjunction  with  the  one-dimensional  transient  heat 
conduction  solution  through  a  spherical  shell  element,  are  ade¬ 
quate  for  predicting  the  temperature  distribution  within  the 
hemispherical  nose  cone.  The  restrictions  are  to  flows  with  a 

Reynolds  Number  range  of  from  one  to  five  million,  stagnation 

0 

temperatures  less  than  4000  R,  and  pressure  distributions  similar 
to  that  given  in  this  report. 
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2.  SYMBOLS 


burled  thermocouple»  located  midway  between  heated  and 
insulated  (inner;  surfaces 

local  skin  friction  coefficient 

Blasius  non-dimensional  friction  factor 

enthalpy 

thermocouple  located  on  heated  surface 
thermocouple  located  on  insulated  surface 
thermal  conductivity 


Reynolds  number*  p  (h  ) 

1/2  *•  ** 


1/2 


"r». 


Reynolds  number  based  on  momentum  thickness* 

Nusselt  number*  Q.  Ra/(T  -  T  )  k 

e  *e 


e/p. 


pressure 
Prandtl  number 

heat  transfer  per  unit  area  per  unit  time 

perpendicular  distance  from  surface  to  centerline  of 
hemisphere 

external  nose  radius  of  hemisphere 


distance  along  surface  measured  from  the  stagnation 
point  along  a  meridian 

temperature 

(T.  -  T)/(I,  -  Ti) 


u 

U 


« 
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e 

p 

p 

T 


local  velocity  at  edge  of  boundary  layer 
uAh, 

e 

pressure  gradient  parameter 
boundary  layer  displacement  thickness 

momentum  thickness 
viscosity  coefficient 

mats  density 

p./p. 

e 

shear  stress 


Subscripts 

aw  adiabatic  wall 

e  conditions  external  to  the  boundary  layer 

i  initial  conditions 

s  stagnation  conditions 

w 


wall  conditions 


4. 


3.  INTRODUCTION 


This  report  is  one  of  a  series  dealing  with  the  general 
problem  of  temperatures «  deformations,  and  material  changes  of 
a  hemispherical  configuration  in  a  hypersonic  flow  environment. 

The  principal  purpose  of  this  report  is  to  present  a  simple  method 
for  predicting  the  temperature  distribution  within  the  hemispheri¬ 
cal  body  under  certain  restricted  hypersonic  flow  conditions.  A 
reasonably  accurate  knowledge  of  the  temperature  distribution  will 
be  necessary  for  subsequent  theoretical  studies  of  stresses  and 
deformations*  and  for  the  interpretation  of  experimental  results. 

First*  a  theoretical  analysis  is  made  of  the  heat  trans¬ 
fer  to  a  hemispherical  body  under  hypersonic  conditions  for  laminar, 
turbulent*  and  transition  boundary  layer  regions.  The  results  are 
compared  with  experimental  heat  transfer  data  obtained  on  a  hemi¬ 
spherical  nose  cone  configuration  in  the  Polytechnic  hypersonic 
facility. 

Except  for  the  region  of  highest  heat  input*  the  agree¬ 
ment  between  theory  and  experiment  is  sufficiently  satisfactory  to 
suggest  the  following  simplification  in  the  temperature  distribu¬ 
tion  estimation  within  the  hemisphere.  The  local  theoretical  heat 
transfer  coefficient  is  used  in  conjunction  with  the  transient 
one-dimensional  heat  conduction  solution  for  spherical  shell  ele¬ 
ments.  Temperature-time  histories  at  selected  points  in  the  in¬ 
terior  are  compared  with  available  test  data  and  the  agreement  is 


5. 


found  to  be  quite  satisfactory  for  purposes  of  future  estimation 
of  temperatures  on  similar  models  which  are  not  equipped  with 
thermocouple  instrumentation. 

Future  hemispherical  nose  cone  models  of  various 
thicknesses*  but  of  similar  external  configuration  will  be  in~ 
strumented  only  on  the  inner*  or  insulated  surface.  Thus*  if 
tne  method  of  estimating  temperatures  presented  in  this  report 
is  reasonably  accurate*  these  predictions  can  be  used  to  obtain 
the  physical  properties  of  materials  from  data  available  on  the 
variation  of  these  propert'as  with  temperature*  to  calculate 
thermal  stresses*  and  for  any  purpose  requiring  a  knowledge  of 
the  temperature  to  make  calculations  and  to  evaluate  experimen¬ 


tal  data. 


} 
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4.  THEORETICAL  HEAT  TRANSFER  CALCULATIONS 

The  heat  transfer  to  a  hemispherical  configuration  un¬ 
der  hypersonic  conditions  was  calculateo  i'cr  tnics  values  of  t '.s 
Reynolds  Number,  «  0.97  x  10^,  3.3  x  10^  and  5.1  x  lO^.  The 

Reynolds  Numbers  were  chosen  to  coincide  with  the  values  at  which 
hypersonic  tunnel  tests  had  previously  been  run  on  a  stainless 
steel  hemisphere.  This  stainless  steel  model  was  specifically  de¬ 
signed  to  obtain  experimental  values  of  the  heat  and  load  input. 
Complete  data  are  reported  in  Ref.  [l],  but  a  brief  description  of 
the  model  and  tests  are  given  in  Section  5. 

The  technique  for  predicting  the  laminar  heat  transfer 
under  the  above  conditions  may  readily  be  obtained  from  Ref.  [2], 
[3],  and  [4],  In  Ref.  [2],  Lees  provides  the  most  convenient  me¬ 
thod  for  estimating  the  laminar  heat  transfer  under  hypersonic  con¬ 
ditions.  It  involves  the  assumption  of  a  constant  value  of  pres¬ 
sure  gradient  parameter,  P,  in  the  transformed  coordinate  system; 
a  ratio  of  wall  to  stagnation  enthalpy,  small  compared  to  unity, 
and  a  linear  dependence  of  viscosity  on  temperature  (p  ja  =p  ^  )* 
There  results  a  simple  numerical  evaluation  of  an  integral  of  ex¬ 
ternal  flow  properties  and  radius  to  determine  the  heat  transfer 
di stribution. 

In  terms  of  the  variables  used  here,  the  theory  of  Lees 


gives 


7 


N  r  -T  puiir 

(0.353)  (»  )  - r 

^  [5  p  u  li  ds] 

0 

The  limiting  process  for  obtaining  the  stagnation  point 
value  from  the  above  equation  leads  to  the  axially  symmetric 
stagnation  point  heat  transfer  (p— ^1,  u*— ♦J’T  etc.)*  This 

gives 

(0.706)  P^^ 

where 


I  =  (du/ds)  at  s  =  0  . 


The  technique  involving  the  determination  of  the  tran* 
sitional  and  turbulent  heat  transfer  requires  the  determination 
of  the  transition  point  on  the  model.  values  of  200  and  300 

were  chosen  for  this  model.  These  transition  Reynolds  Numbers 
are  based  on  experimental  results  performed  by  various  invest!* 
gators . 

Before  starting  the  analysis  of  the  turbulent  boundary 
layert  it  is  necessary  to  evaluate  the  laminar  coefficient  of 
wall  friction  at  the  transitional  point. 


By  definition,  the  wall  friction  coefficient  is: 


Through  Lees  transformation, 


Cf/2  = 


u,  R 

o 


It 

fw  • 


which  may  be  written  as 


Cj/2  =  0.22/N^ 


e 


To  evaluate  the  wall  friction  coefficient  in  the  tran¬ 
sitional  and  fully  turbulent  flow  regimes,  the  friction  law  de¬ 
veloped  in  Ref.  [5]  was  adopted. 

_  0.013  B 

--  J 


N, 


"0  "0 

This  law  is  similar  in  form  to  that  developed  by  Persh 
Ref.  [6].  The  value  of  n  was  taken  as  1,  as  suggested  in  Ref. 
[9]. 

The  constant  B  of  this  law  is  evaluated  from  the  lami¬ 
nar  flow  such  that  the  friction  coefficient  0^/2  remains  continu 
ous  at  the  laminar  to  turbulent  transition  point. 


0.013 
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trans . 
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In  order  to  obtain  the  heat  transfer*  a  stepwise  solu> 
tlon  of  the  boundary  layer  momentum  equation  is  obtained  through 
the  use  of  boundary  layer  form  factors  which  are  based  on  a  cor* 
relation  of  incompressible  data. 

d®  (fc*/®)  +  2  du  1  dp  1  dr 

—  =  (C  /2)  -  9[ - ^  + - + - ] 

ds  u  ds  p.  ds  r  ds 

e  '^e 

In  this  work*  following  the  suggestion  of  Ref.  [5]*  the 
form  factor  was  taken  ast 

“e  =  -  -  1  . 

The  heat  transfer  is  calculated  from 

Q.  =  -  K.)  P, -.(C,/!) 

and  finally 

n  •  n  9 
*e  **  *e 

The  theoretical  distribution  of  heat  transfer  along  a 
meridian  for  a  laminar  boundary  layer  is  shown  in  Fig.  1.  The 
heat  transfer  parameter  correlates  the  heat  transfer  for 

any  value  of  the  Reynolds  Number.  In  the  transition  and  turbu¬ 
lent  boundary  layer  range*  however*  the  calculations  were  made 
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for  three  specific  values  of  the  Reynolds  Number.  Fig.  2  presents 
the  results  of  calculations  based  on  a  transitlon-turbulent  boun- 
dary  layer  for  Reynolds  Numbers  of  0.97  x  10^»  3.3  x  lo^  and 

5.1  X  10^.  Transition  from  laminar  boundary  layer  is  assumed  to 
occur  at  a  Reynolds  Number  based  on  a  momentum  thickness  (N-  )  of 

r  V5  ® 

250.  The  heat  transfer  parameter  N  /N^  is  plotted  against  dis~ 

U  K 

tance  along  the  meridian  from  the  stagnation  point.  In  Fig.  3» 


the  same  calculatlonsy  based  on  a  transition  value  of 
are  graphically  presented. 
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5.  COMPARISON  WITH  EXPERIMENTS 


11  . 


In  Ref.  [l]t  the  results  of  6  test  runs  at  simulated 
hypersonic  speeds  on  a  hemispherical  model  were  reported  in  both 
graphical  and  tabular  form.  The  model  was  fabricated  of  Type 
304  stainless  steel,  had  a  7-3/4  inch  outside  diameter  and  one- 
inch  wall.  It  was  extensively  instrumented  with  pressure  taps 
and  thermocouples.  The  principal  mission  of  these  tests  was  to 
obtain  basic  data  on  both  the  aerodynamic  and  thermal  loads. 
Aerodynamic  load  data  were  presented  in  Ref.  [l]  and  are  also 
shown  in  Fig.  4.  (This  pressure  distribution  was  used  in  cal¬ 
culating  the  theoretical  heat  transfer).  The  procedure  for  ob¬ 
taining  the  thermal  input,  based  on  the  thermocouple  temperature¬ 
time  histories,  is  outlined  below. 

The  hemispherical  model  was  Instrumented  with  a  total 
of  40  thermocouples,  as  shown  in  Fig.  5.  Thirteen  of  these 
thermocouples  were  specially  constructed  one-dimensional  plugs 
for  heat  transfer  work.  This  instrument  was  first  described  in 
Ref.  [7],  and  further  discussed  in  Ref.  [l].  A  sketch  of  a 
typical  thermocouple  plug  is  shown  in  Fig.  6.  The  principal  vir¬ 
tue  of  these  plugs  is  the  relative  ease  with  which  the  heat  in¬ 
put  could  be  deduced  from  the  temperature-time  history  of  the 
thermocouple  on  the  heated  surface  of  the  model.  The  remaining 
thermocouple  installations  were  "standard",  and  were  disposed  on 
the  outer  and  inner  surfaces  of  the  model  as  well  as  at  interior 
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points.  Skotchos  of  a  "standard**  heated  surface  thermocouple  and 
a  burled  thermocouple  Installation  are  shown  In  Figs.  7  and  8. 

The  experimental  values  of  the  heat  transfer  were  de¬ 
duced  from  the  one-dlmenslonal  plug  temperature-time  hlstorlesi 
and  from  similar  data  from  the  standard  thermocouple  Installations 
on  the  outer  (heated)  aurface  of  the  model.  In  the  caae  of  the 
one-dlmenslonal  plugs*  the  heat  conduction  equations  for  the 
transient  one-dlmenslonal  flow  of  heat  through  a  finite  slab  were 
solved  for  the  heat  input  q^  on  the  heated  side  of  the  slab  as  a 
function  of  time.  The  other  aide  of  the  slab  was  considered  to 
be  Insulated.  The  plug  temperature-time  hlatory  was  used  as  a 
boundary  condition.  Except  for  the  first  few  seconds*  the  values 
of  q^  thus  calculated*  resulted  In  an  essentially  constant  value 
of  the  Nusaelt  Number.  The  time  average  of  the  almost  constant 
value  of  the  Nuaselt  Number  waa  used  as  the  experimental  value 
presented  In  this  report. 

In  obtaining  the  experimental  heat  transfer  from  the 
standard  heated  surface  thermocouples*  use  was  made  of  the  tran¬ 
sient  heat  conduction  solution  for  the  flow  of  heat  through  a 
spherical  shell  element*  one-dimenslon.al *  outside  surface  heated* 
Inside  surface  Insulated*  Ref.  [8].  For  the  geometry  parameter 
corresponding  to  the  tested  model*  temperature-time  plots  were 
carefully  drawn  for  constant  values  of  the  heat  transfer  coeffi¬ 
cient.  Nhen  the  experimental  temperstura-tlme  history  was  super¬ 
imposed  upon  this  set  of  curves*  the  heat  transfer  coefficient 
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could  oatily  bo  intorpolatod  at  aavaral  pointa  and  tha  raaulta 
avaragad.  Againi  thaaa  haat  tranafar  coafficiant  valuaa  wara  aa** 
aantially  conatant  aftar  tha  firat  faw  aaconda  of  tha  taat  run. 

A  comparlaon  of  the  axper Imental  valuaa  of  tha  haat 
tranafar  thua  daducad  with  tha  thaoratical  valuaa  la  ahown  in  Figa. 
9f  10  and  11. 

In  Fig.  9»  for  axampla*  tha  thaoratical  curvaa  ara  glvan 
In  the  form  of  Nuaaalt  Number  varaua  dlatanca  along  tha  marldlan 
from  tha  atagnatlon  point  for  a  value  of  N|^  =  0.97  x  10”.  Thaaa 
curvaa  ara  aaally  obtained  from  tha  basic  data  praaantad  in  Figa.  1» 

2  and  3.  Note*  from  Table  1»  that  a  value  of  N|^  =  0.97  x  10”  cor** 
raaponda  cloaaly  to  taat  runa  1  and  2»  and  hence  axparlmantal  data 
from  thaaa  runs  ara  ahown.  Tha  axparlmantal  data  daducad  from  tha 
ona*dimansional  pluga  ara  consistently  higher  than  those  data  da* 
ducad  from  tha  standard  surface  thermocouple  installations.  Near 
tha  stagnation  pointt  there  is  a  discrepancy  of  up  to  2b9(  between 
theory  and  exparimentf  but  there  is  an  improvement  with  higher  v.a- 
lues  uf  s.  Unf ortunatelyt  it  is  not  altogether  certain  whether 
transition  occurs  at  the  lowest  Reynolds  Number>  0.97  x  10^  (Fig.9)t 
but  if  transition  does  occur  it  would  be  at  a  value  of  7  >  0.9. 

Fig.  10  shows  the  comparison  between  theory  and  experi* 

ment  for  =  3.3  x  10^,  corresponding  to  conditions  of  test  runs 

~  6 

3  and  4.  A  similar  presentation  for  =  5.1  x  10  »  corresponding 
to  test  run  5  is  shown  in  Fig.  11.  Once  again  it  can  be  noted 
that  experimental  heat  transfer  values  from  tha  thermocouple  plugs 
are  usually  higher  than  the  heat  transfer  values  daducad  from  the 
standard  surface  thermocouple  data.  In  Figs.  10  and  11  there  is 
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no  question  of  transition  of  the  boundary  layer  from  laminar  to 
turbulent.  Additional  data  in  the  region  of  highest  heat  trans** 
far  would  have  been  desirable.  Several  of  the  instruments  in 
this  areat  howeveri  were  damaged  during  the  installation  of  the 
model  into  the  wind  tunnel  and  during  the  course  of  the  test 
runs . 

On  the  basis  of  the  comparisons  of  Figs.  9»  10  and  11« 
it  was  decided  to  take  the  heat  transfer  results  calculated  for 
the  laminar  case  and  for  N.  -  300  and  use  the  heat  conduction 
equations  through  a  spherical  element  to  predict  the  temperatures 
within  the  hemispherical  model. 


6.  TEMPERATURE  DISTRIBUTION 
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Tht  ttaptratur*  dlttrlbutlon  within  th«  hemitphtr ical 
modal  was  calculatad  by  utilizing  the  haat  conduction  aolutiona 
for  tha  tranaiant  ona-dimanaional  flow  of  haat  in  a  spharical  shall 

As 

alamonti  Rtf.  [8].  For  a  givon  a  and  N^f  th#  haat  tranafar  coaffl- 
ciant  was  takan  from  Figs.  9*  10»  or  11.  Tht  innar  surfaca  was 
aasumad  to  ba  inaulatad.  Tamparatura-tima  historias  could  thara* 
fora  ba  obtainad  at  any  point  within  tha  modal*  for  any  of  tha  thraa 
Raynolds  Numbar  rangas  for  which  tha  thaoratical  haat  transfar  cal¬ 
culations  wars  mada. 

In  tha  six  tast  runs  of  Raf.  CO*  tharnocouplas  wars  bu- 
riad  at  various  points  within  tha  body  (saa  Fig.  5)  and  on  tha 
innar  surfaca.  Tharafora*  tha  points  salactad  for  computing  tha 
tamparatura'Hima  historias  by  tha  abova  mathod*  includad  points  at 
which  tharmocouplas  wars  physically  located  on  tha  modal.  Tha  ob¬ 
vious  purpose  was  to  obtain  a  comparison  batwaan  tha  axparimantal 
and  tha  calculatad  tamparatura-tima  historias. 

Tha  results  of  this  comparison  are  prasantsd  in  Figs.  12 
and  13»  corresponding  to  tha  lowest  and  highest  Reynolds  Numbar* 
raspactivaly .  An  examination  of  these  tamparatura-tima  historias 
reveals  that  tha  calculatad  tamparaturas  are  usually  higher  than 
those  measured  during  tha  tast.  (T  is  defined  in  such  a  manner 
that  it  will  decrease  with  increasing  tamparatura. )  Maximum  devi¬ 
ation  from  tha  axparimantal ly  determined  tamparaturas  occurs  at 
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s  =  1.4|  6.35li  dlfftrence  at  Nr  =  0.97  x  10^,  and  13.4?^  difference 

~  6 

at  Nr  =  5.1  X  10  .  At  the  points  of  highest  heat  transfer*  the 

percentage  error  was  somewhat  less.  For  example,  for  runs  1  and  2, 

~  6  — 

Npj  =  0.97  X  10  ,  thermocouples  24  and  36  at  s  =  0.698  gave  tempera- 
ture>time  histories  for  which  the  maximum  deviation  from  the  calcu¬ 
lated  temperature-time  history  was  approximately  2.5^.  At  the 
highest  Reynolds  Number,  the  comparable  error  at  the  same  location 
was  10^.  It  may  be  noted  that  while  the  maximum  percentage  error 
occurs  at  the  region  of  highest  theoretical  heat  input,  the  maximum 
percentage  errors  in  temperature  occurred  elsewhere. 

The  agreement  between  the  experimental  and  calculated 
temperature-time  histories  is  good  when  one  considers  that  the  cal¬ 
culated  temperatures  are  based  on  a  one-dimensional  heat  conduction 
solution.  The  fact  that  the  calculated  values  of  the  temperature 
are  almost  always  higher  than  the  experimental  values  is  not  wholly 
unexpected  inasmuch  as  the  calculated  heat  transfer  results  are  gen¬ 
erally  higher  than  those  deduced  from  the  plugs  and  surface  thermo¬ 
couples. 
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7.  CONCLUSIONS 

On  the  basis  of  the  comparison  between  theory  and  ex¬ 
periment  of  heat  transfer  and  temperature  distribution  presented 
In  this  report,  it  is  concluded  that  the  simplified  method  used  to 
predict  temperatures  within  a  hemispherical  nose  cone  in  a  hyper¬ 
sonic  flow  is  satisfactory.  The  restriction  must  be  made  to  flows 
in  the  Reynolds  Number  range  of  from  1  x  10  to  5  x  10  ,  and  to 
stagnation  temperatures  less  than  4000®R .  In  order  to  apply  the 
method  it  is  first  necessary  to  calculate  the  theoretical  heat 
transfer  for  the  appropriate  Reynolds  Number.  Using  the  local  val 
ues  of  the  heat  transfer,  In  conjunction  with  the  transient,  one¬ 
dimensional  heat  conduction  solution  through  a  spherical  shell  ele 
ment,  yields  the  temperature-time  history  at  any  point  wlthtM  I'.r 
body. 
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